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ABSTRACT 

This paper presents multi-band photometric follow-up observations of the 
Neptune-mass transiting planet GJ 436b, consisting of 5 new ground-based tran- 
sit light curves obtained in May 2007. Together with one already published light 
curve we have at hand a total of 6 light curves, spanning 29 days. The analysis 
of the data yields an orbital period P = 2.64386 ± 0.00003 days, mid-transit 
time [HJD] = 2454235.8355 ± 0.0001, planet mass Mp = 23.1 ± 0.9 = 
0.073 ± 0.003 Mjup, planet radius Rp = 4.2 ± 0.2 i?e = 0.37 ± 0.01 Rjup and 
stellar radius Rg = 0.45 ± 0.02 Rq. Our typical precision for the mid transit 
timing for each transit is about 30 seconds. We searched the data for a possi- 
ble signature of a second planet in the system through transit timing variations 
(TTV) and variation of the impact parameter. The analysis could not rule out a 
small, of the order of a minute, TTV and a long-term modulation of the impact 
parameter, of the order of -1-0.2 year~^. 
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Introduction 



Of the almost 300 extrasolar planets discovered to datqj, the ~ 35 transiting planets 
are the only ones allowing a measurement of their mass and radius and the study of thei r 



atmospheres fe.g.. ICharbonneau et all 120071 . lFortnevll2008l . lGuillotll2008l . IPont et al.ll2008al ). 



Of those, GJ 436b serves as a unique opportunity to study a planet with mass and radius 
as small as Neptune. Although several planets have already been dete cted with minimum 



masses similar to the mass o f GJ 436b, for ex ample GJ 581 b. c fc d. (jBonfils et al.l 12005 



Udrv et al.ll2007h . HD 4308b (lUdrv et al.ll2006h and Cnc 55e (iFischer et al.ll2008h . GJ 436b 
is the only transiting Neptune-mass planet discovered so far. Moreover, it is the only known 
transiting planet orbiting an M-type star, presenting an opportunity to help constrain stellar 
models in this mass range. 



GJ 436b was initially di scovered bv iButler et al.l (120041 ) through a radial velocity (RV) 
modulation of its host star. Maness et al. (j2007l ) have refined its orbital elements, specif- 
ically the eccentricity, e = 0.160 ± 0.019, and identified a linear velocity trend of 1.36 ± 
0.4 m s~^ yr~^- Those authors suggested that the velocity trend might result from the pres- 
ence of a long-period second planet in the system, a planet that could induce a periodic 
modulation of the orbital eccentricity of GJ 436b. This typ e of effect was suggested already 



i n a b inary system for 16 Cygni B (IMazeh et al. 



for triple stellar systems ( Mazeh fc Shaham Il979, see also Mazeh 1990), a n d for a planet 



1997; Holman et al. 



19971 ). iManess et al. 



(120071 ) pointed out that this interpretation of the eccentricity looks especially attractive if 
the circularization timescale of GJ 436b is shorter than the age of the system, because then 
we need to explain why the orbit of the planet has not been circularized. 



Butler et al.l (120041 ) obtained photometric observations of GJ 436 at the expected time 



of possible transits and concluded that a transit w as unlikely. Neverth eless, the transit- 
ing nature of GJ 436b was recently discovered by iGillon et al.l (l2007al ) who measured a 
planetary radius of ~ 4 i?® and mass of 22.6 M®. Soon after the discovery, a transit and 



secondary echpse events were observed by Sp itzer Space Telescope at 8 //m (IGillon et al. 



2007bl : iDeming et al.ll2007l : iDemory et al.ll2007l ). further constraining the system parameters 
and reducing the uncertainties by a factor of about 10 for the planet to star radius ratio and 
mid-transit timi ng. Those obs ervations resulted in a somewhat increased planetary radius. 



of 4.2 ± 0.2 iTorred (120071 ) was able to refine further the planet parameters by deriving 



more accurate constraints on the host star and obtained a slightly larger planetary mass, 
23.17 ± 0.79 M®, and a more precise planetary radius, of 4.22^ 



)+0.09 p 
-0.10 



^for an updated list see: |http: / /exoplanet.eu/ 
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Recently, iRibas et al.l (120081 ) restudied the system and suggested that the observed radial 
velocities of the system are consistent with an additional small, super- Earth planet in a close 
orbit around G J 43 6 in a 1:2 mean-motion resonance with the known planet. Similarly to 
Maness et al.l (120071) . they suggested that this planet could pump eccentricity into the orbit. 
Ribas et al. ( 20081 ) further suggested that such an additional planet can induce a precession 
of the orbital plane of GJ 436b, if the orbital planes of motion of the two planets are inclined 
relative to each other. Such a precession should change the inclination relative to our line 
of sight and therefore might explain why iButler et al.l (120041 ) failed to observe a transit in 
2004. 



Immediately after the detection of the transits of G J 436b by iGillon et al.l (l2007af ) we 
launched a ground-based (GB) observational campaign to obtain high-quality transit light 
curves in different filters. Our motivation was to obtain the first group of light curves which 
will be used by future studies to look for a possible variatiori of the impact parameter and to 
search for transit timing variations (TTV; lAgol et al.ll2005l : iHolman fc Murrayl l2005l ) . Our 
observations are described in §[2j In our analysis, described in §[3], we simultaneously analyze 
six complet e GB transit l ight cu rves: five light curves obtained in our campaign and the light 
curve from lGillon et al.l (j2007al ). In § |l]we discuss our results. 



2. Observations 



In addition to the transit of UT 2007 May 2, which was reported by lGillon et al.l (j2007al ) 
and was the first complete transit observation of GJ 436b, we obtained five complete transit 
light curves of four different transit events. The following paragraphs briefly describe these 
observations. 

We observed the transits of UT 2007 May 4 and May 25 with the Wise Observatory 1 m 
telescope, located in Israel. We used a Princeton Instruments (PI) VersArray camera and a 
1340 X 1300 pixels back-illuminated CCD, with a pixel scale of 0'.'58 per pixel, giving a field 
of view of 13'0 x 12'6. On May 4 we used a Johnson V filter and on May 25 a Bessel R filter. 
The flexible scheduling of the Wise 1 m telescope was an important factor in obtaining the 
transit light curve of May 4, less than three days after the first observation of a complete 
transit. 

The transit of UT 2007 May 4 was simultaneously observed with the Wise Observatory 
0.46 m telescope, operated remotely. The camera was a 2148 x 1472 pixels Santa Barbara 
Instrument Group (SBIG) ST-10 XME CCD detector with a field of view of 40^5 x 27^3 
and a scale of I'.'l per pixel. This camera has no filters. For a detailed description of this 
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telescope and instrument see iBrosch et al.l (120081 ) . 



On UT 2007 May 23 and May 31 we used the 1.2 m telescope at the Fred L. Whipple 
Observatory (FLWO) on Mount Hopkins, Arizona. The camera was the KeplerCam, which 



is a 4 096^ Fairchild 486 CCD with a square field of view 23' 1 on a side (jSzentgyorgyi et al. 



20051 ). On both nights the filter was Sloan z and a 2 x 2 binning was used, giving a pixel 



scale of 0'.'68 per binned pixel. 

Basic data reduction procedures, including bias and fiat-field correction, were applied to 
all images using standard IRAF@ routines. We then performed aperture photometry of GJ 
436 and several comparison stars of similar brightness showing no significant variability. The 
light curve of GJ 436 was calibrated by dividing it by the summed fiux of the comparison 
stars. Next, we fitted a polynomial of degree one or two to the out-of-transit (OOT) points 
vs. time, and divided all points by this polynomial. The amplitude of these polynomial 
corrections in mmag per transit duration was in the range of 0.1-1.0 mmag. As a final 
step, we divided the entire light curve by the OOT median intensity and normalized the 
measurement uncertainties so the median uncertainty OOT will equal the OOT RMS. Tabled] 
lists all photometric measurements of the 5 light curves obtained in this campaign. 

Table [2] lists for each light curve the UT start date and time, observatory and telescope, 
the filter used, limb darkening coefficients used in its analysis, average exposure time, average 
cadence (in min~^), duration of the entire observation, start and end airmass, the /5 factor 
(see below), RMS residuals, and the Photometric Noise Rate (PNR). The PNR is a quantity 
which takes into account the RMS and the cadence, defined as RMS x \J cadence, in units of 



mmag per minute (see also lBurke et al.ll2008l . section 2.2). 



3. Data Analysis and Results 



We decided not to include the publicly available Spitzer 8 /xm light curve in our simulta- 
neous fitt ing since the shape of that light curve is dependent on how the Spitzer "ramp" was 
modeled (IGillon et al.ll2007bl : iDeming et al.l 120071 ). As the Spitzer measurements are much 
more precise than the GB measurements, this ramp could have systematically infiuenced our 
results. For comparison, the Spitzer light curve PNR is a factor of ~3 smaller than that of 
our ground-based light curves. 



^The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 
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Accounting for correlated noise (jPont et al.l |2006| ) was done similarly to the "time- 
averaging" method of lWinn et al. J2008f i7 After a preliminary analysis we binned the residual 
light curves using bin sizes close to the duration of ingress and egress. The presence of 
correlated noise in the data was quantified as the ratio between the binned residual light 
curves standard deviation and the expected standard deviation in the absence of correlated 
noise. For each light curve we estimated (3 as the largest ratio among the bin sizes we used. 
We then multiplied the relative flux errors by f3. Values of /3 are in the range of 1.0-2.1 and 
are listed in Table [H 



3.1. Simultaneous fitting of all parameters 



We fitted the system parameters to all six GB light curves simultaneously. Our transit 
light curve model consisted of nine parameters: The orbital period, P; a particular mid- 
transit time, Tc] planet to star radius ratio, Rpl Rf,\ semimajor ax is scaled by the stellar 
radius, a/Rg] impact parameter, b (See Eq. 3 of IWinn et al.l 120071 for the formula of the 
impact parameter in an eccentric orbit); two limb-darkening coefficients ui and U2, for a 
quadratic limb-darkening law; orbital eccentricity, e; and longitude of periastron, u. The 
latter four parameters were held fixed in the fitting process, hence our model had 5 free 
parameters. 

For transiting planets on eccentric orbits the mid point, in time, between transit start 
and end is different than the time of sky-projected star-planet closest approach. For GJ 436b 
this difference is comparable to the typical uncertainty on the mid-transit times obtained 
here. In our model the mid-transit time is defined as the time of closest approach, which is 
also the time of minimum light. 



For e and u we used the values given bv lManess et al.l (120071 ). Although other authors 



(e.g., iDemory et al.l 120071 . iDeming et al.l 120071 ) report more precise values for e, their uncer- 
tainty on uj is either much higher than that of iManess et al.l (120071 ) or not reported. As the 
errors in two orbital parameters are correlated, we chose to adopt the reference giving the 
smallest uncertainty on both. To detenmnejii and U2 we adopted a model of Teg = 35 00 K 
jManess et al1l2007n Jg^= 4.84 3 JTorresl l2007h and [Fe/H]=-0.03 JSonfils et al.lboosh and 
interpolated on the I Clare tl (120001 . |200J) grids. For the light curve obtained with the Wise 
Observatory 0.46 m telescope, which has no filter, we used limb-darkening co efficients cor 



respo nding to the R filter, as its CCD response resembles a "wide-i?" filter (IBrosch et al. 
20081 ). The coefficients used for each light curve are listed in Table [2l 



We used the formulas of iMandel &: Agoll (120021 ) to determine the relative fiux as a func- 
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tion of the plane t-star sky-proiected s e paration and the Monte Ca rlo Markov Chain (MCMC) 



algorithm (e.g.. iTegmark et al.l 12004 



FordI I2OO5I : iGregoryl |2005| ) to determine the parame- 



ters that best fit the data, along with their uncertainties. MCMC algorithm s have already 



been used extensively in the literature for fitting transit light curves (e.g., iHolman et al. 
20061 : ICoUier Cameron et al.l I2OO7I : Icillon et all boosi : ISurke et all I2OO7I ) . Assuming a uni- 
form prior for the five fitted parameters, the algorithm can be viewed as a random walk in 
the five-dimensional parameter space of the model. Starting from an initial guess, at each 
step the algorithm examines a new point in the parameter space which is reached by adding a 
Gaussian random value to each of the parameters. The algorithm decides whether to accept 
the new point in the parameter space and move there, or repeat the current point in the 
chain depending on the resulting posterior probability. If the new point has a higher relative 
posterior probability the algorithm will move there, and if not, it will move to the new point 
with a relative probability of exp(— Ax^/2), where Ax^ is the difference between the two 
points in parameter space. The distribution of parameter values in all points is used to 
derive the best fit value and its uncertainties. Here we took the distribution median to be 
the best fit value and the values at the 84.13 and 15.87 percentiles to be the +la and —la 
confidence uncertainties, respectively. 

The widths of the Gaussian distributions used to find a new point are a crucial part of 
the MCMC alg orithm and they affect the fraction of accepted points, /. In order to have 
/ close to 25% ( IGregoryl I2OO5I : iHolman et al.ll2006l ) we divided each long chain, of 500,000 
steps, into smaller chains (mini-chains), of 1,000 steps, and re-evaluated the distribution 
widths after the execution of each mini-chain. This re-evaluation was done by scaling them 
according to the relative difference between / of the previous mini-chain, and the target 
value of 25%. The same scaling was applied to all parameters. 

To estimate the relative size of the distribution widths for each one of the model param- 
eters we initially performed a sequence of 10 mini- chains , where only that parameter was 
allowed to vary and the target / was 50% (IGregoryl |2005| ). Thus, a single MCMC run was 



comprised of several sequences of 10 mini-chains and a long chain, of 500 mini-chains. For 
determining the final result we considered only the long chain while ignoring the first 20% 
of its steps. 

We performed 10 MCMC runs as described above while changing each time the initial 
guess and initial distribution widths. Results of those runs were highly consistent with each 
other and our final result is based on all runs. 

Before describing our results we briefiy describe how we tested our analysis: 



To test the adaptive form of our MCMC algorithm we re-ran it while fixing the distri- 
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bution widths during the long chain part. The differences in the resuhs for the fitted 
parameters did not exceed O.lScr, and the differences in the estimated uncertainties 
were no more than 10%. 

We repeated our analysis while skipping the polynomial corrections mentioned in Sec- 
tion [21 The results for the fitted parameters were less than 0.2a from the original ones, 
and the estimated uncertainties were similar at the 20% level. 



We applie d our analys i s on th e Spitzer 8 /ini 



similar to iGillon et al.l fl2007bf ). lOeming et all (120071 ) and ISouthworthI fl2008f ). at the 



ight curve alone, and d erived results 



~ lo" level. 

To further test the validity of our results we re-ran MCMC while adopting a Gaussian 
prior pro bability for e and u . The central values and standard deviations were those 
given by lManess et al.l (120071 ). The results were not modified significantly. 



Fig. [T] presents the six GB light curves, with our fitted model overplotted in solid lines. 
The three light curves obtained in the beginning of 2007 May are of better quality, especially 
the Euler 1.2 m V and Wise 1 m y light curves, with RMS residuals of 1.2 and 0.9 mmag, 
respectively. The increased scatter in the three light curves obtained at the end of 2007 May is 
the result of the high airmass of these observations, done close to the end of the observational 
season. 

Our estimates for the fitted parameters, while fitting all six light curves simultaneously, 
are listed in Table [31 For comparison, this table includes results from previous studies when- 
ever these parameters were fitted directly. In case they were not, but can be calculated from 
other fitted parameters, they are given without an error. iGillon et al.l (l2007al ) based their pa- 
rameter determinations on an analysis of the Euler 1.2 mV light curve from UT 2007 May 2, 



which is included in this work. The other th ree studies referred to in Table [31 (IGillon et al. 



2007bl : [Peming et al 



12007 



Southworthll2008l ) are based on the Spitzer 8 /xm light curve. Re- 



sults obtained by [Torres et al.l (120081 ) are not presented since for the light curve parameters 
of GJ 436b these authors gave the weighted average of previous works. 

Table [H lists our results for the physical parameters of the system, along with the 
corresponding val ues obtained previously. Assuming the mass of t he star to be — 



0.452 ± O.O13M0 ( lTorresll2007l ) and adopting the orbital parameters of iManess et al.l (120071 ) 
(specifically: K, e and u) we derived the planet mass and semi-major axis. Using Newton's 
version for Kepler's Third Law we derive the star and planet radius. 
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3.2. Fitting transit timing and impact parameter for each light curve 



As mentioned in § [T|, iRibas et al.l (120081 ) suggested the presence of a second planet or- 
biting GJ 436, in an outer orbit close to a 2:1 mean- motion resonance with the transiting 
planet. The gravitational interaction between th e two planets may c ause a detected TTV 
signal (IHolman fc Murravl l2005l : lAgol et al.ll2005l ). iRibas et al.l (120081 ) suggest that the sec- 
ond planet is responsible for changing the inclination angle, of the transiting planet orbit. 
The presence of such a planet may be detected through measuring transit timing and the 
impact parameter of many transit light curves. As the light curves obtained here are the first 
group of such light curves we performed two additional MCMC analyses, described below. 

First, we fitted Rp/Rg, a/Rg and h to all light curves and an independent Tc to each 
transit light curve, in order to estimate the mid-transit time of each event. In this analysis 
the period does not determine the time intervals between transit events. It is used only 
for determining the true anomaly of each measurement within each light curve so it is only 
weakly constrained. Hence, in this analysis P was fixed at the value derived in our original 
analysis. The resulting six values were later used to derive the transit ephemeris by a 
linear fit. The new period is ~ O.SScr from the period fitted in our original MCMC run. 
Repeating the process described in this paragraph with the new period results in a variation 
of only 0.02(7, so no further iterations were done. 

Table [5] contains the mid-transit timings fitted independently to each light curve. Note 
that light curves with the smallest RMS residuals do not necessarily have the smallest mid- 
transit time uncertainty. For example, the FLWO 1.2 m light curve, at = 3, has relatively 
high RMS residuals although a small mid-transit time uncertainty. This is simply the effect 
of the short KeplerCam cadence, as the increased number of points act to better constrain 
the mid-transit time. This effect is quantified by the PNR (see Table [2]). 



By fitting a linear function. 



T,{E) = T,(0) + P X 



to the transit epoch, E, we derived P = 2.64386 ± 0.00003 days and T,(0) = 2454235.8355 ± 
0.0001 [HJD]. The transit event with E = was chosen to be right in the middle of our 
observed events, although we did not observe that particular event. We take this ephemeris 
as our final result, listed in Table [3l and note that it is consistent with P and Tc derived 
from our original MCMC fitting. 

The residuals from the linear fit are presented in Fig. known as the 0-C diagram. 
The figure shows some indication for a variability of the period, refiected by the high 
value of the fit, of 23 for 4 degrees of freedom. However, the small number of points does not 
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allow to assign any significance to the detection of this variability, especially because some 
systematic effects could shift some of the points. 

Second, we performed a MCMC analysis where we looked for a transit impact parameter 
variation. The four parameters P, T^, Rp/Rs and a/Rg were fitted to the entire data while 
the impact parameter, b, was fitted independently to each light curve. The derived impact 
parameters are listed in Table [5] and plotted in Fig. [3l As a reference, we over-plotted in 
Fig. [3] the value of our original run, in a solid line, and the upper and lower la confidence 
limits in dashed lines. No significant variation can be seen in the data, although a long-term 
trend of the order of 0.2 year~^ can not be ruled out. 

The two separate analyses described in this section were repeated without applying the 
polynomial corrections mentioned in Section [21 The results were consistent with the original 
ones. In addition, for both the 0-C transit timing residuals and the impact parameters 
derived here we looked for a possible correlation with several external parameters, including 
mid-transit airmass, filter central wavelength and the light curves polynomial correction 
amplitude. All correlations were below the 1.5a confidence level. 



4. Discussion and Summary 



This paper presents an analysis of 6 complete transit light curves of GJ 436b, spanning 
29 days, or 11 orbital periods. 



Our period is consistent with the period given by iManess et al.l ( 120071 ) . of 2.64385 ± 
0.00009 days, and is more precise. Our value for Rp/Rg (see Tabl e ED, is large r, by 1- 2 
g, than the r a dii ra tio d erived for the Spi tzer 8 fim light curve by iGillon et al.l (l2007bl ). 
Deming et al.l (120071 ) and ISouthworthI (120081 ). If this difference is real it may be induced by 
stellar spots, or wavelength depende nt opacities in the planetary atmosphere (lBarmanll2007l . 
Pont et al.ll2008al . iPont et al.ll2008bl ). Interestingly, a similar effect, i.e., an increased radius 



ratio in the optical relative to the IR, was already observed for HD 189733b (jPont et al. 
2007h . 



Our derived physical parameters, mainly the planetary radius and mass, are consistent 
with those derived previously (See Table S]). 

We examined the individual transit timing of each light curve (TTV, see Fig. [2]) and 
derived the impact parameter of each light curve independently (see Fig. [3]). The transit 
timing residuals show a possible hint for variability, supported by the high of the fit- 
However, many more light curves are needed to establish the variability and to explore its 
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nature. The impact parameter plot shows a possible trend, although a linear fit to the six 
data points gives a slope consistent with zero. 



Comparing our de rived transit ephemeris to transit times published recently flAlonso et al 



20081 . iBean et al.ll2008l ) shows that the predicted transit times are somewhat earlier than the 
measured times. The difference is at t he 1-2(7 lev e l, so i t is not highly significant. However, 
if real, it supports the claim made by iBean et al.l (l2008l ) for a long-term drift in the transit 
times. Together with the long-term slope in GJ 436 RVs (IManess et al.l 120071 1 it suggests 
the existence of another object in the system, with a long period. Future photometric and 
spectroscopic data will be able to better study this possibility. 

As mentioned in the introduction, the impact parameter of GJ 436b can be changed by 
a precession of the line of nodes, induc ed by a second planet whose orbital plane is inclined 



relative to t he orbit of GJ 436 b (e. g., Miralda-Escude 20o3). Such an effect was already 



suggested by lSoderhjelml (119751 ) and lMazeh fc ShahamI (Il976i ) for stellar triple systems. The 
precession period is 

^ 1 P (2) 



node 



P,2 



( iMiralda-Escudel |2002| ). where Mp2 and P'?, are the inass an d period of the second planet, 
respectively. For the planet suggested by iRibas et al.l (|2008| ) , this precession is of the order 
of 10 years. 

For eccentric orbits the impact parameter can also be modified by the precession of the 
periastron, b ecause the impac t parameter depends on the orientation of the line of apsides 
(see Eq. 3 of I Winn et al.l 120071 ). The apsidal motion can be driven by the stellar quadrupole 
moment, or by a second planet. The precession driven by the stellar quadrupole moment 
could be quite slow, and therefore probably does not contribute to the variation of the impact 
parameter. The rate of the apsidal precession induced by another planet is of the same order 
of magnitude as the rate of the nodal precession, estimated in Eq. [2l and therefore can be 
of the order of years, depending on the parameters of the unseen additional planet in the 
system. 

In order to estimate the possible implication of the apsidal precession on the impact 
parameter we plot in Fig.Hlthe impact parameter of GJ 436b as a function of the argument of 
the periastron tu, at fixed orbital eccentricity and inclination. We can see that h is modulated 
between ~ 1.0 and 0.7 over the apsidal precession period. The upper dashed line of the figure 
indicates that at some phase of the precession period the transit of GJ 436b will even lose 
its fiat-bottom shape. The figure suggests that if the apsidal motion is of the order of ten 
years, we should be able to observe soon a change of the impact parameter. In case this 
modulation is due to apsidal motion, we will then be able to confirm the change of the line 
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of apsides direction by timing the secondary eclipse or by precise RV measurements. 
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the Kepler Mission under NASA Cooperative Agreement NCC 2-1390 with the Smithsonian 
Astrophysical Observatory. 



REFERENCES 

Agol, E., Steffen, J., Sari, R., & Clarkson, W. 2005, MNRAS, 359, 567 

Alonso, R., Barbieri, M., Rabus, M., Deeg, H. J., Belmonte, J. A., & Almenara, J. M. 2008, 
A&A, 487, L5 

Barman, T. 2007, ApJ, 661, L191 

Barnes, J. W. 2007, PASP, 119, 986 

Bean, J. L., et al. 2008, A&A, 486, 1039 

Bonfils, X., et al. 2005, A&A, 443, L15 

Bonfils, X., Delfosse, X., Udry, S., Santos, N. C, Forveille, T., & Segransan, D. 2005, A&A, 
442, 635 

Brosch, N., Polishook, D., Shporer, A., Kaspi, S., Berwald, A., & Manulis, I. 2008, Ap&SS, 
314, 163 

Burke, C. J., et al. 2007, ApJ, 671, 2115 
Burke, C. J., et al. 2008, ApJ, 686, 1331 

Butler, R. P., Vogt, S. S., Marcy, G. W., Fischer, D. A., Wright, J. T., Henry, G. W., 
Laughlin, G., & Lissauer, J. J. 2004, ApJ, 617, 580 

Charbonneau, D., Brown, T. M., Burrows, A., & Laughlin, G. 2007, in Protostars and 
Planets V, B. Reipurth, D. Jewett, & K. Keil (eds). University of Arizona Press 
(Tucson), p. 701 (|arXiv:astro-ph/0603376[) 



- 12 - 



Claret, A. 2000, A&A, 363, 1081 
Claret, A. 2004, A&A, 428, 1001 
Collier Cameron, A., et al. 2007, MNRAS, 380, 1230 

Deming, D., Harrington, J., Laughlin, G., Seager, S., Navarro, S. B., Bowman, W. C, & 
Horning, K. 2007, ApJL, 667, L199 

Demory, B.-O., et al. 2007, A&A, 475, 1125 

Fischer, D. A., et al. 2008, ApJ, 675, 790 

Ford, E. B. 2005, AJ, 129, 1706 

Ford, E. B., Quinn, S. N., & Veras, D. 2008, ApJ, 678, 1407 
Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661 
Fortney, J. J. 2008, ArXiv e-prints, larXiv:0801.4943| 
Gillon, M., et al. 2007, A&A, 472, L13 
Gillon, M., et al. 2007, A&A, 471, L51 

Gillon, M., Triaud, A. H. M. J., Mayor, M., Queloz, D., Udry, S., & North, P. 2008, A&A, 
485, 871 

Gregory, P. C. 2005, ApJ, 631, 1198 

Guillot, T. 2008, Physica Scripta Volume T, 130, 014023 

Holman, M., Touma, J., & Tremaine, S. 1997, Nature, 386, 254 

Holman, M. J., & Murray, N. W. 2005, Science, 307, 1288 

Holman, M. J., et al. 2006, ApJ, 652, 1715 

Mandel, K., & Agol, E. 2002, ApJ, 580, L171 

Maness, H. L., Marcy, G. W., Ford, E. B., Hauschildt, P. H., Shreve, A. T., Basri, G. B., 
Butler, R. P., & Vogt, S. S. 2007, PASP, 119, 90 

Mazeh, T., Krymolowski, Y., & Rosenfeld, G. 1997, ApJ, 477, L103 

Mazeh, T. 1990, AJ, 99, 675 



- 13 - 

Mazeh, T., & Shaham, J. 1976, ApJ, 205, L147 

Mazeh, T., & Shaham, J. 1979, A&A, 77, 145 

Miralda-Escude, J. 2002, ApJ, 564, 1019 

Pont, F., Zucker, S., & Queloz, D. 2006, MNRAS, 373, 231 

Pont, F., et ah 2007, A&A, 476, 1347 

Pont, F., Knutson, H., Gilhland, R. L., Moutou, C, & Charbonneau, D. 2008a, MNRAS, 
385, 109 

Pont, F., Gilhland, R. L., Knutson, H., Holman, M., & Charbonneau, D. 2008b, to appear 
in MNRAS f arXiv:0810.5731) 

Ribas, I., Font-Ribera, A., & Beaulieu, J.-P. 2008, ApJ, 677, L59 

Seager, S., Kuchner, M., Hier-Majumder, C. A., & Militzer, B. 2007, ApJ, 669, 1279 

Soderhjelm, S. 1975, A&A, 42, 229 

Southworth, J. 2008, MNRAS, 386, 1644 

Szentgyorgyi, A. H., et al. 2005, Bulletin of the American Astronomical Society, 37, 1339 
Tegmark, M., et al. 2004, Phys. Rev. D, 69, 103501 
Torres, G. 2007, ApJ, 671, L65 

Torres, G., Winn, J. N., & Holman, M. J. 2008, ApJ, 677, 1324 
Udry, S., et al. 2006, A&A, 447, 361 
Udry, S., et al. 2007, A&A, 469, L43 
Winn, J. N., et al. 2007, ApJ, 665, L167 
Winn, J. N., et al. 2008, ApJ, 683, 1076 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



-14- 



Table 1: Photometry of GJ 436b. The table includes the 5 light curves obtained in this work 
and will appear in entirety in the electronic version of the journal. 



Observatory+ 


Filter 


HJD 


Relative 


Relative 


Telescope 






flux 


flux error 


Wise 0.46m 


Clear 


2454225.229864 


1.0015 


0.0020 


Wise Im 


V 


2454225.218116 


0.9997 


0.0010 


FLWO 1.2m 


z 


2454243.658856 


1.0002 


0.0020 


Wise Im 


R 


2454246.309414 


1.0024 


0.0017 


FLWO 1.2m 


z 


2454251.641074 


0.9987 


0.0024 



Table 2. List of light curves analyzed in this work. 



Start Date 




Start Time 


Observatory+ 


Filter 


Ul 


U2 


exp. 


Cadence 


Duration 


start 


end 


13 


RMS 


PNR=' 




UT 




[hh:mm] UT 


Telescope 








time [s] 


[min^-'] 


[hours] 


AM 


AM 




[mmag] 


[mmag min^-'-] 


2007 


May 


2 


00:08 


Euler 1.2 m 


V 


0.340 


0.444 


80 


0.44 


3.86 


2.10 


2.24 


1.0 


1.2 


1.8 


2007 


May 


4 


17:39 


Wise 0.46 m 


clem' 


0.343 


0.398 


20 


1.23 


3.19 


1.03 


1.14 


1.6 


2.0 


1.8 


2007 


May 


4 


17:22 


Wise 1 m 


V 


0.340 


0.444 


60 


0.65 


2.90 


1.05 


1.07 


1.5 


0.9 


1.1 


2007 


May 


23 


03:56 


FLWO 1.2 m 


z 


0.088 


0.522 


10 


2.44 


3.91 


1.02 


1.99 


1.3 


2.0 


1.3 


2007 


May 


25 


19:33 


Wise 1 m 


R 


0.343 


0.398 


40 


0.88 


3.60 


1.15 


3.79 


1.3 


2.4 


2.6 


2007 


May 


31 


03:30 


FLWO 1.2 m 


z 


0.088 


0.522 


10 


2.38 


3.20 


1.02 


1.64 


2.1 


2.1 


1.4 



^Photometric Noise Rate = RMS x VCadence. 



I 

I— ' 

I 



Table 3: Light curve fitted parameters. 



Reference 


P 


Te-2454200 


Rp/Rs 


a/Rs 


h 




[d] 


[HJD] 








Tliis work 


2.64386 


35.8355 


0.085 


13.6 


0.85 




±0.00003 


±0.0001 


±0.001 


±0.5 


±0.01 


Gillon et al. 2007ai 






0.082 












±0.005 






Gillon et al. 2007b^ 




80.78148 


0.0830 




0.849 






+0.00015 






+0.010 






-0.00008 






-0.013 


Deming et al. 2007^ 




80.78149 


0.0839 


13.2 


0.85 






±0.00016 


±0.0005 


±0.6 


+0.03 












-0.02 


Southworth 2008^ 




80.78174 


0.08284 


13.68 








±0.00011 


±0.00090 


±0.51 





1 Period was fixed at the Maness et al. 2007 value, of P = 2.64385 ± 0.00009 days. 



Table 4: Results for the physical parameters. 



Reference 


a 






Rs 




Rp 




[AU] 


[Me] 






This work 


0.02872 +0.00030 


23.1 ±0.9 


0.45 


±0.02 


4.2 


±0.2 




-0.00025 












Gillon et al. 2007a 




22.6 ±1.9 


0.44 


±0.04 


3.95 


+0.41 














-0.28 


Gillon et al. 2007b 






0.463 


+0.022 


4.19 


+0.21 










-0.017 




-0.16 


Deming et al. 2007 


0.0291 ±0.0004 


22.24 ±0.95 


0.47 


±0.02 


4.33 


±0.18 


Torres 2007 


0.02872 ±0.00027 


23.17 ±0.79 


0.464 


+0.009 


4.22 


+0.09 










-0.011 




-0.10 


Southworth 2008^ 






0.452 


±0.017 


4.08 


±0.16 



^Assuming = 0.4521°;°^^ Mq (Torres 2007), orbital parameters of Maness et al. (2007) 
and the period derived here. 
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Table 5: Mid-transit time and impact parameter fitted independently to each light curve*. 



E 


Te - 2454200 




t 




b 


Ab 


t 




[HJD] 


[min] 












-5 


22.61612±0.00037 


-0.15 ■ 


-0.28 


0.; 


838±0.017 


-0.013 


-0.61 


-4 


25.26002±0.00038 


-0.10 ■ 


-0.18 


0.: 


828±0.020 


-0.023 


-0.96 


-4 


25.26052±0.00030 


0.62 


1.45 


0.: 


843±0.016 


-0.008 


-0.37 


3 


43.76657±0.00026 


-0.71 ■ 


-1.88 


0.: 


846±0.018 


-0.005 


-0.23 


4 


46.40982±0.00040 


-1.58 ■ 


-2.74 


0.; 


843±0.014 


-0.008 


-0.41 


6 


51.69956±0.00030 


1.34 


3.12 


0.; 


854±0.016 


0.003 


0.13 



* The table gives the result of two separate MCMC analyses. In each analysis only a single 
parameter (Tc or b) was fitted separately to every light curve and all other parameters 
were fitted to all light curves simultaneously, except for P which was fixed when fitting an 
individual to each light curve. 

^ Difference between the measured Tc and the expected mid-transit time according to the 
fitted ephemeris (see Table [3]). The difference is given in minutes and also in units of the 
uncertainty on each mid-transit time. 

* Difference between the measured impact parameter while fitting each light curve separately, 
and the result of the simultaneous fit (see Table ED- The difference is given also in units of 
the uncertainty on each impact parameter. 
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1.01 



>< 1 



0.99 

a 
<u 

0.98 



2007 May2-Euler1.2m V 



1.01 

1 

0.99 
0.98 



2007 May 4 - Wise 0.46m clear 



22.5 22.55 22.6 22.65 



25.25 25.3 25.35 



1.01 



>< 1 



0.99 

DC 

0.98 



2007 May 4 - Wise 1 m V 



1.01 
1 

0.99 
0.98 



2007 lyiay 23 - FLWO 1 .2m. z,^^. 




25.2 25.25 25.3 25.35 



43.7 43.75 



43.8 



1.01 



>< 1 



0.99 

m 
o 
CC 

0.98 



2007 May 25 - Wise 1 m R 



' . ■/'"■•■•A •■.•V". 



46.3 46.35 46.4 46.45 
HJD - 2454200 



1.01 
1 

0.99 
0.98 



2007 May 31 - FLWO 1 .gm..? 




51.65 51.7 51.75 
HJD - 2454200 



Fig. 1. — The six light curves analyzed in this work. The overplotted solid line is our best 
estimate model. For each light curve, residuals from the model are plotted, centered on a 
relative flux of 0.98 for clarity. Within each panel the UT date, observatory, telescope and 
filter used are given. 
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Fig. 2. — Observed minus calculated (0-C) transit timing of the six light curves included in 
this study. The graph shows the residuals from a linear fit to the transit timing as a function 
of the transit epoch. Table [5] lists the actual transit timings. 
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0.87 



0.86 



0.85 



>0.84 



C) 



o 



o 



() 



o 



0.83 



() 



0.82 



0.81 



0.8 





Epoch 



Fig. 3. — The impact parameter, b, fitted independently to each hght curve, as a function 
of the transit epoch. The sohd hne marks the best estimate derived in our original MCMC 
run, while fitting the impact parameter to all light curves simultaneously. The dashed lines 
mark the upper and lower la confidence limits. Table [5] lists the actual values of the impact 
parameters presented here. 
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1.1 




0.65 - 



0.6' ' ' ' ' ' ' ^ 

50 100 150 200 250 300 350 

CO [deg] 

Fig. 4. — The solid line presents the impact parameter, b, vs. u, for the eccentricity and 
inclination of GJ 436b. The filled circle (green) marks the current position, error bars are 
too small to be shown. The upper and lower dashed lines (red) mark the impact parameter 
when the stellar radius is multiplied by a factor of 1 — Rp/Rg and 1 + Rp/Rg, respectively. 
When multiplying the stellar radius by 1 — Rp/Rg (upper dashed line) the impact parameter 
will be smaller than one only for non-grazing transits. When multiplying by 1 + Rp/Rg 
(lower dashed line) then the impact parameter will equal one when the minimal star-planet 
sky projected distance equals the radii sum, i.e., for 6 > 1 there will be no transit at all. 



22 



